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Summary. We present the current status of our construction of synthetic spectra 
for type la supernovae. These properly take into account the effects of NLTE and an 
adequate representation of line blocking and blanketing. The models are based on a 
sophisticated atomic database. We show that the synthetic spectrum reproduces the 
observed spectrum of 'normal' SN-Ia near maximum light from the UV to the near- 
IR. However, further improvements are necessary before truly quantitative analyses 
of observed SN-Ia spectra can be performed. In particular, the inner boundary 
condition has to be fundamentally modified. This is due to the dominance of electron 
scattering over true absorption processes coupled with the flat density structure in 
these objects. 



1 Introduction 

A detailed understanding of the physics and the explosion mechanism of type la 
supernovae (SNe la) is essential to minimizing systematic errors of cosmological 
parameters derived from the interpretation of SNe la observations. A key role in 
this context play realistic models of the expanding atmospheres of SNe la which 
are based on the full non-equilibrium physics involved, which link the results of 
current hydrodynamic explosion models (see, e.g., Reinecke et al. 2002 [9], Gamezo 
et al. 2003 [1]) to the comprehensive observational data. 

The spectra of SNe la contain useful information about the energetics of the 
explosion — luminosity and velocity — and the nucleosynthesis, especially in shells 
above the pseudo-photosphere, which are primarily observable in early epochs. To 
gain access to this information several groups have started to compute synthetic 
optical spectra for SNe la in the photospheric epoch (Pauldrach et al. 1996 [8], 
Nugent et al. 1997 [5], Lentz et al. 2001 [3], Hoflich et al. 2002 [2]); so far, however, 
all use simplifying assumptions with uncertain influence on the resulting spectrum. 
Realistic models for quantitative analysis are still missing. 

The main reason we still do not have realistic models is that the physical con- 
ditions in expanding SN la-atmospheres — high radiation energy density in a low 
matter density environment — make this modeling especially difficult. Simplifying 
assumptions such as local thermodynamic equilibrium (LTE) are invalid because ra- 
diative processes dominate over local collisional processes. Thus, a realistic model 
requires at minimum a solution of the full non-LTE problem. Other difficulties 
arise because the ejecta do not contain significant amounts of hydrogen or helium, 
elements that form the main contribution to continuous opacity in other stellar 



2 D.N. Sauer, A.W.A. Pauldrach, T. Hoffmann, and W. Hillebrandt 



objects. This causes a strong dominance of line opacity over continuous opacity, 
which is further complicated by the high expansion velocities of the ejecta (up to 
~ 30000 km/s). Thus, thousands of Doppler-shifted and -broadened spectral lines 
of low-ionized metal atoms overlap and form a non-thermal "pseudo-continuum" 
with strong line-blocking and line-blanketing effects. 

Furthermore, the entire light emission is powered by 7-photons (and in later 
times also positrons) originating from the decay of 56 Ni — > 56 Co — > 56 Fe. These 
7-photons are redistributed to lower energies by scattering and ionization events. 
They are unlikely to be thermalized — in best case only partially — even within the 
pseudo-photosphere and the outer-most parts of the atmosphere because the density 
increases only moderately toward the center. Recent explosion models predict that 
56 Ni can be present even at high velocities (Reinecke et al. 2002 [9]). This indicates 
the need for a consistent treatment of the energy input within the atmosphere, which 
may significantly influence the emergent spectrum. Also, the observed changes in 
the spectra with time as the photosphere recedes into the ejecta show that the 
abundance of the elements varies strongly with radius. 

All these factors have to be taken into consideration to develop a realistic ra- 
diative transfer model that allows a reliable quantitative analysis of SNe la. 

In this paper, we describe the computation of synthetic spectra for SNe la 
based on a consistent treatment of the full NLTE radiative transfer (see Sauer & 
Pauldrach 2002 [10]). This method will be used to test hydrodynamic explosion 
models against the observations and will eventually provide a powerful tool for 
spectral diagnostics of SNe la. 



2 The Method 

The basic approach for modeling supernova spectra is analogous to the treatment 
of O-stars with expanding atmospheres because the physics involved is similar in 
many respects. Our code development is based on a program package which has 
proven to be very successful for the quantitative spectral analysis of hot stars with 
radiation-driven winds (Pauldrach et al. 2001 [7], Pauldrach 2003 [6]). 

The models treat the expanding atmosphere as stationary and spherically sym- 
metric. The underlying density structure, velocity field, and composition of the 
ejecta (averaged over angles) are provided by the explosion model. 

The code provides a consistent solution of the radiative transfer and the NLTE 
rate equations for the occupation numbers of the atomic levels. A proper treatment 
of line blocking and blanketing effects is included. All significant contributions to 
opacity and emissivity from bound-bound, bound-free, and free-free transitions, as 
well as Thomson-scattering, are taken into account — for this purpose we have com- 
piled and computed a comprehensive atomic database involving accurate radiative 
and collisional ionization and excitation cross-sections and atomic level energies of 
all important ions. 

After each iteration cycle, the temperature for each radius point is determined 
by balancing all energy gains and losses to the gas — implementation of the 7-photon 
deposition is currently under development; thus far, it is only considered with the 
assumption that all deposition takes place below the photosphere. Together with 
the use of the diffusion approximation at the the lower boundary this implies the 
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assumption of a complete thermalization of the 7-photons. The total luminosity 
(given by the amount of synthesized 56 Ni and the 7-deposition efficiency) is a pa- 
rameter that is used to fit the shape of an observed spectrum. 

Finally, the emergent spectrum is given by the computed radiation field at the 
outer boundary of the converged model. 



3 First Results 

Fig. 1 (see Sauer & Pauldrach 2002 [10]) compares the synthetic spectrum of a 
NLTE-model based on the W7 explosion model (Nomoto et al. 1984 [4]) to the 
observed spectrum of SN 1992A. At first glance the spectra seem to agree quite 
well. However, a closer look reveals a few discrepancies: the model reproduces the 
UV and blue part of the observed spectrum fairly well, but it fails to match the 
relatively steep decline toward infrared bands. Additionally, the absorption features 
of the model that are deeper than the observed may suggest that the peak fluxes 
are too large. 




Fig. 1. Synthetic spectrum for a SN la around maximum light based on the W7 
explosion model with homogenized composition. For comparison, the observed spec- 
trum of SN1992A is also plotted. 



Thus, we cannot yet consider the models to be sufficient for quantitative spec- 
troscopy of SNe la. Relevant SN la physics, in particular the treatment of the 
7-deposition and a better analytical approach to the boundary condition at the 
photosphere (see next section) still have to be included into the model. 
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4 Discussion and Conclusions 

The enhanced fluxes of the model redward of ~ 4000 A shown in Fig. 1 result from 
the invalid assumption of a photosphere analogous to the atmospheres of hot stars. 
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Fig. 2. The contibutions to the opacity in a stellar atmosphere (upper panel) and a 
supernova (lower panel). Shown are the in both cases the Thomson-opacity XThom 
(the grey straight line) and the true continuum opacity Xcont (the dark solid line). 
Xtruo denotes the sum of the continuum opacity and the line opacity. 

Various conditions in the ejecta contribute to the breakdown of the diffusion- 
approximation at the lower boundary. For stellar atmospheres, the diffusion limit for 
radiation is an excellent approximation as the conditions deep in the photosphere 
are very close to LTE. This is because the opacities at all wavelengths are dominated 
by true processes: H and He provide a strong bound-free continuum and free-free 
processes are supported by the exponential density increase at the photosphere (see 
Fig. 2). 

In SNe la, however, the opacity is dominated by electron scattering over a large 
wavelength range, between the optical red and the radio band (see Fig. 2) — the 
sources of true continuous opacity (free-free and bound-free) are orders of magni- 
tude smaller because the composition does not contain hydrogen and helium, and 
the density distribution is flat. The electron scattering opacity, in contrast to true 
opacities, does not couple the radiation field to the local temperature. As the flux at 
the inner boundary from the diffusion approximation is determined by the tempera- 
ture gradient which is constrained only by the true processes and does not consider 
scattering dominance. The input flux at these wavelengths is overestimated. We 
note these conditions do not change at larger optical depths. Moreover, the sit- 



Synthetic Spectra for Type la Supernova 



5 



uation already occurs at early epochs where the assumption of a photosphere is 
generally considered to be a valid approximation. 

These results show that the commonly used lower boundary condition (i.e., 
analytical expressions based on the diffusion-approximation) is highly problematic 
for SNe la. An inner boundary that correctly accounts for the physical conditions 
is required in order to correctly solve for the radiation field in SNe la. 

Overall, the first results obtained with our method are encouraging. The main 
spectral features of SN la spectra are well reproduced, indicating that the underly- 
ing physics is correctly described. However, the results also show that even at early 
epochs the SNe la envelopes are comparable neither to the atmospheres of hot stars 
nor to gaseous nebula. Thus, the assumptions usually made for these cases cannot 
be applied to the expanding envelopes of SNe la. In particular, a valid analytic 
approach to the boundary condition at the inner radius of a scattering-dominated 
atmosphere has to be formulated. 
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